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Abstract 



We discuss an elementary problem in electrostatics: What does the charge dis- 
tribution look like for a free charge on a strictly one-dimensional wire of finite 
length? To the best of our knowledge this question has so far not been discussed 
anywhere. One notices that a solution of this problem is not as simple as it might 
appear at first sight. 
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I Introduction 



Some time ago one of the authors of this paper was confronted with the following 
question by a student: What does the charge distribution look like for a free charge 
on a strictly one-dimensional wire of finite length, if we consider the usual Coulomb 
repulsion law? Despite its rather trivial appearance this problem led to controversial 
discussions. The immediate "obvious" answer by most considering the problem for 
the first time is to suggest a charge distribution like the one in Fig. 1: That is one 
intuitively expects an accumulation of charge at the ends of the wire since there is no 
repelling charge outside. We urge the reader to make up his mind too before reading 
on. 

Let us remark that one-dimensional systems of electrons interacting with long- 
range Coulomb forces have attracted much attention recently since it has become pos- 
sible to realize them experimentally as one-dimensional semiconductor structures 3 . 
To avoid misunderstandings we want to emphasize that this paper does not aim at 
contributing in this direction of research, though this serves as an interesting back- 
ground. Quantum effects and the fermionic nature of electrons are not taken into 
account here. We rather want to introduce a nice exercise in classical electrostatics 
for students or anybody else interested in elementary problems. Only very elementary 
tools of mathematics and physics will be used, still at first sight the answer might 
appear counter-intuitive and surprising. 

Let us make one more remark in order to avoid confusion. We do not discuss a 
problem in a one-dimensional world, but rather a one-dimensional problem defined 
in three-dimensional space. Therefore we use the three-dimensional Coulomb law 
V(r) oc |r| _1 . Notice that starting from the one-dimensional Laplace equation the 
interaction potential has the different form V(r) oc \r\. This then defines the well- 
investigated one-dimensional Coulomb gas model that has been solved exactly in the 
thermodynamic limit independently by Prager 4 and Lenard 5 . Unfortunately the 
techniques used in these solutions cannot be carried over to our problem 6 and we have 
to rely on other tools in the sequel. 

II Regular izat ion prescriptions 

Once one starts analyzing the problem, one immediately notices that it is ill-defined 
in its original formulation. For simplicity we will assume that the wire has unit length 
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throughout this work: The charge distribution p(x) is defined on the interval [0, 1] of 
the x-axis. In order to solve our problem we want to minimize the energy functional 



E \p] = o / dx i d V 
I Jo Jo 



\x - y\ 

under the constraints of a fixed total charge Q here set to one 



/ dx p(x) = 1 
Jo 



and 

p(x) > 0. (3) 
The integral (|l|) diverges since 

X ~ E + f 1 ) dy , 1 , ^° 21ne. (4) 
o Jx+eJ \x — y\ 

We will discuss two obvious possibilities to make the problem well-defined, i.e. ways 

to regularize the energy functional (H). It is not immediately clear that both give the 

same answer. 

(i) Define the one-dimensional wire to be the limit of ellipsoids in three-dimensional 
space when two semiprincipal axes of the ellipsoid shrink to zero. This way we 
use a regularization by going to a well-defined problem in the embedding three- 
dimensional space, thereby avoiding the singular one-dimensional problem. This 
approach is discussed in Sec. III. 

(ii) A physically appealing approach is to put n equally charged classical particles 
in equilibrium on the wire. Each particle has an individual charge = 1/n. 
This defines the looked for charge density in the limit n — > oo. Although this 
regularization procedure requires more effort than (i), it is more convincing in 
the sense of being a "microscopic" approach. This makes up the main part of 
our paper and is worked out in Sec. IV. 

Of course other regularization procedures are also possible. A very natural choice 
would be 

Vd(r) = -™L= ±3 V(r) = ± (5) 
Vr 2 + d 2 \r\ 

since this is usually used for one-dimensional systems interacting with long-range 

Coulomb forces 7 . The main advantage is that the one-dimensional Fourier transform 

of Vd(r) exists for finite d. However, no exact solution seems possible for finite d and 

it is therefore difficult to investigate the limit d — > 0. 
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Ill Shrinking ellipsoids 



Consider the conducting ellipsoid in Fig. 2 with semiprincipal axes a, b and c. One 
can argue that for b, c — > the ellipsoid shrinks to a one-dimensional wire of finite 
length 2a. 

Now the potential problem of a free charge Q on a conducting ellipsoid has been 
well-known for a long time and is treated in many advanced textbooks on electrostatics, 
see for example Ref. . The explicit solution relies on the fact that the Laplace equation 
is separable in elliptic coordinates. Let us simply quote the result for the potential 



v(0 



Q 



dX 



to J Z + a?) (A + b 2 ) (A + c 2 ) 

with £ = £(x, y, z) defined implicitly by 



(6) 



x 



+ 



+ 



i. 



(7) 



a? + £ b 2 + £ c 2 + £ 
The surface charge a is given by the normal derivative of the potential at the surface 



a = — 



dn 



5=o 



In rectangular coordinates this is 
Q 1 



a 



47r abc 



I x 2 y 2 z 2 



(9) 



For our purposes we can assume b = c. One would argue that the total surface charge 
in the strip S in Fig. 2 collapses onto the line charge t(x) dx of the wire at that point x 

t(x) = <t(x) ■ S(x), (10) 

where S(x) dx is the surface of the strip S. It is a simple exercise in geometry to show 



S(x) dx = 2nb 



\ 



1 



1 



dx. 



On the other hand, evaluating Eq. (0) at the surface gives 
Q 1 



a(x) 



An ab x 2 / 



(12) 
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This implies that independent of b the line charge density defined like in Eq. ([H]) is 
constant along the wire 




According to this reasoning there are no finite size effects in our one-dimensional wire 
problem: There is no accumulation of charge at the ends of the wire! 8 

The sceptical reader can certainly question the validity of this proof by pointing 
out that even for infinitesimal b, in some respects the ellipsoid is in no way similar to 
a wire. E.g. the ends of the ellipsoid are always much thinner than its middle, and the 
smaller available space may compensate a charge accumulation at the ends of the wire. 
Therefore we will use a different, more microscopic regularization in the next section. 

IV Regularization by discretization 

We regularize the singular one-dimensional problem by considering n equally charged 
classical particles with individual charges 1/n in the interval [0, 1]. Due to Coulomb 
repulsion there is a unique equilibrium state with charges at positions In the limit 
n — ■> oo this will define the continuum charge distribution on the wire that we are 
interested in. 

IV. 1 Upper and lower bounds 

First of all we will show that the energy of the discretized problem diverges like Inn 
but can still be determined within bounds of width |. This will be done by calculating 
an upper and a lower bound for the energy of the equilibrium distribution. 

An upper bound is given by the energy E max of a uniform distribution with charges 
at positions Xj = (i — l)/(n — 1), i = 1 . . . n 



We use the well-known asymptotic behaviour (see e.g. formula 0.131 in Ref. 2 ) 




(14) 



n 



1 




t=l 



(15) 
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where 7 = 0.5772 ... is Euler's constant, and find 



n — 1 



n-l 



E max (n) = — — \2^2lnj + 2n^ + 0(lnn) 
n \ 3=1 

= Inn + 7 - 1 + 0(l/n)lnn. (16) 

In the last step we have used Stirling's formula. 

In order to derive a lower bound E m i n we first of all notice that the energy of a 
discrete distribution can be split up as 

E = E 1 + E 2 + . . . + (17) 

where Ek is the sum of energies of all particles to the k-th neighbours (see Fig. 3) 

in 1 

E k = 4 E 1 r ( 18 ) 

i+fc <n 

It is straightforward to see that each is bounded from below by the energy obtained 
if a uniform distribution of the k-th neighbours is assumed. This gives e.g. in Eq. (|18|) 

£i>A>-l) 2 . (19) 
Furthermore we have for even n 

E2>-J^-l] 2 (20) 



n 2 V2 
and for odd n 

1 ( fn + l \ 2 /n-l \ 2 \ 2 /n \ 2 

- 1 ) + (— " x ) J 5 ^ (2 " 1 ) • (21) 

These results are easily generalized for arbitrary fc 



k fn 



^nAl- 1 )- (22 » 

A lower energy bound -E m i n is therefore 

min = ^SHjfc-V = SU"n + ^J 

= lnn + 7- - + 0(l/n) (23) 

where we have again used Eq. (|T5]). 

Comparing Eq. ( P3[ ) with the upper bound from Eq. flTE| ) one notices that we have 
obtained rather strict limits for the energy of the equilibrium state with minimum 
energy. Obviously a uniform distribution cannot be too far from the true equilibrium 
solution with minimum energy. That this is indeed true will be shown in the next 
section. 
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IV. 2 Rigorous results for the charge distribution 

Theorem Suppose the charge distribution p n {x) forn particles converges for n — > oo 
towards a continuous distribution p(x) defined on ]0, 1[. Then p{x) is constant. 

Proof: Let us assume that the theorem is wrong and there exist yx,y2 £ ]0, 1[ with 
p(yx) > p{y-i)i see Fig. 4. It shall be shown that there exists an N so that for n > N 
transferring a particle from yx to the middle of two particles at y 2 results in a decrease 
of the energy. Hence the particles were not in equilibrium for n > N and the theorem 
is proven. 

This can be seen as follows: First of all we regard the energy needed to trans- 
fer a particle P from yx to y 2 while neglecting the particles in [yi — e, yx + e] and in 
[y 2 — £, 2/2 + £]■ Here e has been chosen so small that both intervals are completely in 
]0, 1[ and do not overlap. The change of energy is 

\t(J^--^L) ^ O(e), (24) 




yi - t\ \y 2 - t 

in [yx - e, j/i + e] and in [y 2 -e,y 2 + e] ^ 
pit) elsewhere. 

It is essential to note that Eq. fl24|) does not contain any term divergent in n, but has 
a finite and well-defined limit for n — > 00. Let us now calculate the energy of the 
particle P with respect to the other particles in the interval [yx — e,yi + e]. We can 
assume the distribution to be locally uniform, deviations from that will only contribute 
in order e 2 . Thus the particles in the interval }yx,Vi + e] are at positions 

£ 

Xi = Vi-\ 1 , « = !■■■«?, (26) 

rag 

where g = e p(yi) is the total charge in ]y l7 yi + e]. This gives 

2 nq 1 nepfei) / n 

nAEx = -E^^ = 2 E — + 0(.) 
" i=i Fi-3/i| i= x 1 

= 2p( yi )\n(p(yx)£n) + O(n ) , (27) 

where the sum is done like in the previous section. Particle P is placed in the middle 
between two particles at y 2 and this costs the extra energy AE 2 

nAE 2 = 2 E P -H + 0{e) 
i=i ? 2 

= 2p(y 2 )ln(p(i/2)^) + O(n ) . (28) 
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The total change of energy is 

n(A£ 7 + AE 1 - AE 2 ) = 2 (p( yi ) - p(y 2 )) Inn + O(n ) . (29) 

For n large the expression is dominated by the first term on the r.h.s., i.e. the charge 
distribution is stable only for p(yi) = p{y 2 )- □ 
Two remarks are to be made: 

• The proof does not say anything about the endpoints of the interval. It can be 
shown that if the charge density at the endpoints is well-defined in the continuum 
limit, it will obey the inequality p(0) < 2p(0.5). 

• The proof can be generalized for all potentials V(x — y) = \x — y\~ a with a > 1. 
Therefore only potentials with longer range forces than the Coulomb potential 
(i.e. potentials with a < 1) can show finite size effects in the wire problem. 



IV. 3 Some estimates for the discrete charge distribution 

After proving that the continuum charge distribution is flat, let us go back to the dis- 
cretized version of the problem. We want to derive some estimates for the distribution 
of the charges close to the ends of the wire. In combination with numerical results 
in the next subsection, this will help us to reconcile the somehow counter-intuitive 
picture of a flat continuum charge distribution. 

We use the same notation as in Sec. IV. 1. The distances between particle i and 
particle i + 1 are denoted by di = x i+ i — x^. Then force equilibrium for the second, 
third etc. particle on the wire means that the following set of equations is fulfilled 

111 1 

~d\ ~ ~d\ + {d 2 + d 3 y + {d 2 + d 3 + d A y + • • ' 
1111 i 

+ 32 = 32 + 73 ■ + 73 ■ , ■ + ■ ■ ■ ( 30 ) 



(rf! + rf 2 ) 2 d\ 4 (d 3 + d 4 ) 2 (d 3 + d A + d 5 ) 2 
1 1111 1 

(di + d 2 + d 3 ) 2 (d 2 + d 3 ) 2 dl ~ d\ (d A + d b ) 2 {d A + d 5 + d 6 ) 2 

The first and the n-th particle are trivially at positions x\ — and x n — 1 and no 

equilibrium condition can be formulated for them. We sum the first / of the above 

equations (/ < n/2) and subtract equal terms on both sides 

l j l m j n—2 I 

^ / / 777, \ 2 ^ / ^ / / 777, \ 2 ^ ^ ^ / / 777, \ 2 (^^) 

m=l 
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For simplicity we assume that n is even. Then obviously the largest distance between 
particles is d max = d n / 2 and we have the following inequality 

l i n / J j n-2 \ 

Et^i^ £- + ' E ^ ■ (32) 

m=l y> , "max \m=l '' t m=l+l " L J 

The first sum is like Eq. (|i"5|). The second sum can be evaluated using the integral 
approximation 

„ n-3/2 

n-2 i r 1 1 1 

I T A « / / dm At = ^— i — r - • (33) 

^ m 2 7 m 2 + 1/2 n-3/2 v ; 

m=Z+l l+1/2 I I 

One then finds 

ETT^>-^(^ + 7 + l4) (34) 



m=l 



\ — ' , 1 ^max \ ^ / 



\fc=l / 

for large values of I, n. For I = n/2 one has d\ < d 2 < . . . < d n /2, therefore 

n/2 1 1 n/2 j j ^2 



/ m X 2 < ^2 H m 2 - ^2 « • ( 35 ) 



m=l / y> , \ r M ,„ : (/ T C 

\fc=l 



We finally get 



n 2 



d\ < d 2 max 2 _ (36) 

In n — In 2 + 7 + - 

Lower bounds can be given too using the first equilibrium condition in Eq. (p0|) 
111 1 

dj ~d% + (d 2 + d 3 y + (d 2 + d 3 + d A y + " ' 



1_ ^ 2 1 1 7T 2 

J2 

l 2 m =l 

6 -d* 
r 2"2 



" dl ^ m 2 " rf 2 6 (37) 



d?>-|d°. I* 



Eq. ( p6|) says that for large n, the distance of the first two particles on the wire is 
arbitrarily smaller than the distance of two particles in the middle of the wire. Numer- 
ical calculations in the next subsection will show that the equality is nearly realized 
in Eq. (|38|). Thus the particles close to the ends of the wire show strong finite size 
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effects. However, these effects vanish in the continuum limit as we have proved in the 
last subsection. This discrepancy is resolved by noticing that the finite size region 
observed in the discretized problem shrinks for n —>■ oo as will be demonstrated in the 
next subsection. 

IV. 4 Numerical results 

In order to gain a better understanding of the problem, we have also used numerical 
methods to find the equilibrium position of the n particles sitting on the line of unit 
length. We have employed the Hybrid Monte-Carlo 9 algorithm that updates the posi- 
tions of the individual particles until they find their equilibrium positions: The molec- 
ular dynamics part of the algorithm moves all particles according to the electrostatic 
forces, which are computed only once for all particles at each step. The Monte-Carlo 
update scheme ensures that the algorithm is exact, i.e. that equilibrium is reached. 
During the simulation, the "temperature" was reduced as the charges moved more and 
more closely to their equilibrium positions. Convergence was verified by checking that 
the ground state energy did not decrease any more for additional sweeps within a given 
numerical precision. In general, convergence was very quick since the starting point of 
the simulation with all particles equidistant is quite "close" to the equilibrium position 
for the reasons explained in Sec. IV. 2. Simulations with up to n = 8193 particles have 
been performed. 

The first interesting quantity investigated numerically is the ground state energy 
E(n) as a function of the number of particles n plotted in Fig. 5. Also drawn are the 
upper and lower bounds derived in Eqs. flTED and fl2"3"|). One notices that the measured 
energies seem to settle on the upper bound with good precision for n — > oo. This 
raises the interesting (and so far unanswered) question whether the upper bound flTE|) 
becomes the asymptotically exact result for the ground state energy. 

Let us also emphasize that the data in Fig. 5 paedagogically demonstrate the danger 
of extrapolating to n — > oo based on simulations for finite n, at least for systems with 
long range forces: Extrapolating on the basis of the data for n < 1000 particles the 
upper bound would eventually be violated (see the curve in Fig. 5)! Only for more 
than 1000 particles the curve for E{n) bends over and the extrapolation respects the 
exact results 10 . 

In Fig. 6 the quotient of the smallest (at one end of the wire) to the largest (in the 
middle of the wire) particle distance is plotted as a function of n. The upper bound 
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from Eq. ( PU] ) is respected, in particular the interparticle distance at the end of the wire 
becomes much smaller than the maximum distance. On the other hand, the ratio of 
the distance of the first to the second particle d\ and the second to the third particle 
approaches a constant nonzero value for n — > oo. This is shown in Fig. 7 together 
with the lower bound derived in Eq. fl38|). One clearly sees the strong finite size effects 
already mentioned in the previous subsection. 

However, when one looks at the distribution of all charges along the wire plotted 
in Fig. 8 for various values of n, one notices that the regions with strong finite size 
effects close to the ends of the wire shrink very slowly with increasing n, probably as 
slow as 1/ Jhx(n), as suggested by Eq. (|36|). One can interpret these numerical results 
as follows: Eventually the continuum charge distribution appears flat, although for 
a given n the first few particles do never approach equidistant positions. But in the 
continuum limit n — > oo these non-equidistant regions eventually vanish as compared 
to the rest of the wire. This scenario combines the analytical and numerical results 
established in this work. 



V Summary 

We have discussed the seemingly trivial, at least conceptually simple problem of a free 
charge on a one-dimensional wire of finite length. The first thing to notice was that 
this problem is ill-defined in its original formulation due to the diverging ground state 
energy. The question how the continuum charge distribution looks like can only be 
answered after introducing some regularization procedure. 

Two particularly intuitive regularizations have been employed in this paper; shrink- 
ing ellipsoids in Sec. Ill and charge discretization in Sec. IV. Both led to the answer 
that the continuum charge distribution is flat, that is there are no finite size effects! 
The interested reader will be able to find other regularization prescriptions that lead 
to the same answer. 

A flat charge distribution will also be found for potentials 

V(x-y) oc \x-y\- a (39) 

with shorter range forces than the Coulomb potential, that is generally for a > 1. 
For exponents a < 1 the problem is well-defined and finite without the need for 
regularization. It is easy to show that then there are finite size effects. Thus the 
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Coulomb law is the limiting case between finite size effects and no finite size effects in 
a strictly one-dimensional wire of finite length, see Fig. 9. 

One interesting problem left open in this respect is the analytic form of the con- 
tinuum charge distribution for exponents a < 1. Another question that we have not 
been able to answer is whether the ground state energy as a function of the number of 
particles n in Sec. IV is really given by the upper bound Eq. (jTED plus subdominant 
corrections 0(lnn/n). The numerical results in Fig. 5 seem to indicate this. 
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Figure captions 

Fig. 1. A frequent first guess for the charge distribution on a finite wire of unit length. 
Some readers might also suggest divergencies at the ends of the wire. 

Fig. 2. A free charge on a conducting ellipsoid. 

Fig. 3. Sum of energies to the k-th neighbours. 

Fig. 4. Proof by contradiction — a nonflat charge distribution cannot be stable in the 
limit n — > oo. 

Fig. 5. Ground state energy E(n) for a system of n particles with individual charges \ jn 
interacting with Coulomb potentials on a wire of unit length. 

Fig. 6. The quotient of smallest to largest particle distance as a function of the number 
of particles n. 

Fig. 7. The quotient of the distances of particle #1 and #2 and of particle #2 and #3 
on the wire as a function of n. 

Fig. 8. Charge distribution on the wire for various numbers of particles n. Xi is the 
position of particle j^i. Nearest neighbor distances plotted as a 

function of the position along the wire. 

Fig. 9. Finite size effects on one-dimensional wires of finite length for potentials 
V(x — y) oc \x — y\~ a . 
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